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Abstract. A notion of implicit difference equation on a Lie groupoid is intro- 
duced and an algorithm for extracting the integrable part (backward or/and 
forward) is formulated. As an application, we prove that discrete Lagrangian 
dynamics on a Lie groupoid G may be described in terms of Lagrangian implicit 
difference equations of the corresponding cotangent groupoid T*G. Other sit- 
uations include finite difference methods for time-dependent linear differential- 
algebraic equations and discrete nonholonomic Lagrangian systems, as parti- 
cular examples. 



Dedicated to Ernesto Lacomba on the occasion of his 65th birthday 
1. Introduction 

Lie algebroids and groupoids have deserved a lot of interest in recent years since 
these concepts generalize the traditional framework of tangent bundles and its dis- 
crete version, cartesian product of manifolds, to more general situations. In par- 
ticular, it is well-known that the geometric description of the Euler-Lagrange equa- 
tions of a mechanical system determined by a Lagrangian function _L, relies on the 
intrinsic geometry of the tangent bundle TQ, the velocity phase space of a con- 
figuration manifold Q. In the case when the Lagrangian L is invariant under the 
action of a Lie group G, the description, in this case, relies on the geometry of 
the quotient space TQ/G and the equations describing the dynamics are called 
Lagrangc-Poincarc equations [4 . In this sense, Weinstein [55] showed that the com- 
mon geometric structure of the Lagrange- Poincare equations is essentially the same 
as the one of the Euler-Lagrange equations, namely that of a Lie algcbroid. In the 
case of Lagrangian systems on the usual tangent bundles, the tangent bundle carries 
a canonical Lie algebroid structure which is given by the usual Lie algebra of vector 
fields on Q. In the case of reduced lagrangian systems, we use the Atiyah algebroid 
TQ/G — > Q/G to describe the evolution equations. 

In [IS], it is described geometrically discrete Lagrangian and Hamiltonian Me- 
chanics on Lie groupoids, in particular, the type of equations analyzed include the 
classical discrete Euler-Lagrange equations, the discrete Euler-Poincare and discrete 
Lagrange-Poincare equations. These results have applications for the construction 
of geometric integrators for continuous Lagrangian systems (reduced or not) [9j [24] . 

This work has been partially supported by MICINN (Spain) Grants MTM2009-13383, 
MTM2010-21186-C02-01 and MTM2009-08166-E, project "Ingenio Mathematica" (i-MATH) No. 
CSD 2006-00032 (Consolidcr-Ingenio 2010), the project of the Canary Islands government SOL- 
SUB200801000238 and the European project IRSES-project "Geomech-246981" . D. Iglesias wishes 
to thank MICINN for a "Ramon y Cajal" research contract. 
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On the other hand, in [27], W.M. Tulczyjew proved that it is possible to 
interpret the ordinary Lagrangian and Hamiltonian dynamics as Lagrangian sub- 
manifolds of suitable special symplectic manifolds. These results were successfully 
extended to general Lie algebroids in [T^] . More generally, in [T3J [T5J [TB] implicit dif- 
ferential equations are considered. More precisely, an implicit differential equation 
on a manifold Q is a submanifold D C TQ. Given such a datum, the problem of 
integrability is discussed, looking for a subset S C D where, for any ueS there ex- 
ists a curve 7 : / — > Q such that 7(0) = v and 7(f) G D for any t € I. An algorithm 
for extracting the integrable part of an implicit differential equation is formulated. 
More precisely, if D C TQ is an implicit differential equation, the following sequence 
{(C k ,D k )} is considered 

D° = D, C° = tq(D), 
D k = D k-i n TC k-i t C k = TQ (D k ), 

where tq : TQ —> Q is the canonical projection. Clearly, D k C D k ~ 1 and C k ~ l C 
C k for any k. Under some smoothness assumptions, it is proved in [15] that if the 
algorithm stabilizes, that is, there exists k such that D k = £)*"" 1 , then D k C D is 
the (possibly empty) integrable part of D. The relation with Lagrangian mechanics 
is the following one. Given a Lagrangian L : TP — > R, one can consider the 
implicit differential equation Sl Q T(T*P) using the image of the differential of L, 
dL(TP) C T*(TP) and the canonical isomorphism A M : T*{TP) -> T{T*P). If 
we apply the integrability algorithm, the first step will give us the solution of the 
Euler-Lagrange equations. Indeed, if {q i ,v l ) are local coordinates on TP then 

. i dL i dL 
, a ,1 f , i dL i dL i i d (dL\ dL 

Taking as starting point the results of discrete mechanics on Lie groupoids [T01H8] 
and a remark on implicit differential equations on Lie algebroids in [11] , in this paper, 
submanifolds of a Lie groupoid G are interpreted as systems of implicit difference 
equations on G. We study the problem of integrability of these implicit difference 
equations, proposing an algorithm which extracts the integrable part (backward 
and/or forward) of this type of systems and analyze the geometric properties of 
them. In our opinion, this approach offers greater conceptual clarity. Discrete La- 
grangian systems (unconstrained and nonholonomic) are studied within the frame- 
work of implicit difference equations. Moreover, for the unconstrained case, the 
resulting implicit difference equation is a Lagrangian submanifold, thus obtaining a 
discrete version of the results in |12) . Other examples include finite difference meth- 
ods for time-dependent linear differential algebraic equations. The results about 
implicit difference equations are useful for systems more general than the studied 
on this paper; in particular, discrete optimal control theories, discrete systems with 
external constraints ... See for instance fTUlUZO] . 
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2. Preliminaries: Lie algebroids and Lie groupoids 

In this Section, we will recall the definition of a Lie groupoid and some generalities 
about them are explained (for more details, see [LI]). 

A groupoid over a set M is a set G together with the following structural maps: 

• A pair of surjective maps a : G — s- M, the source, and j3 : G — s- M, the 
target. These maps define the set of composablc pairs 

G 2 = {(g,h)eGxG/l3(g) = a(h)}. 

• A multiplication m : G2 — » G, to be denoted simply by m(g, h) = gh, 
such that 

- a(gh) = a{g) and (3(gh) = (3(h). 

- g(hk) = (gh)k. 

• An identity section e : M — > G such that 

- ^(a{g))g = 9 and ge(/3(g)) = g. 

• An inversion map i : G — > G, to be denoted simply by i(g) = g~ x , such 
that 

- 9~ X 9 = e(P(g)) and 99^ = e(a(ff))- 

A groupoid G over a set M will be denoted simply by the symbol G =4 M. 

The groupoid G z4 M is said to be a £ie groupoid if G and M are manifolds and 
all the structural maps arc diffcrentiable with a and (3 differentiable submersions. 
If G =t M is a Lie groupoid then m is a submersion, e is an injective immersion and 
i is a diffeomorphism. Moreover, if x € M, a~ 1 (x) (resp., (3~ 1 (x)) will be said the 
a-fiber (resp., the (3-fiber) of x. 

On the other hand, if g £ G then the left-translation by g £ G and the right- 
translation by g are the diffcomorphisms 

l g : a -\(5( g )) —> a -\a{g)) ; h — ► Z ff (/i) = gh, 
r g :p-\a{g))^p-\P{g)) ; ft — > r fl (ft) - %. 

Note that Z^ 1 = and r^ 1 = r g -i. 

A vector field X on G is said to be left-invariant (resp., right-invariant) if 
it is tangent to the fibers of a (resp., f3) and X(gh) = (TiJ g )(Xh) (resp., X(gh) = 
(T g r h )(X(g)), for (g,h)£G 2 . 

Now, we will recall the definition of the Lie algebroid associated with G. 

We consider the vector bundle t : AG — > M , whose fiber at a point x £ M 
is A X G = V e ( x )a = Ker(T e / x \a). It is easy to prove that there exists a bijection 
between the space of sections of r, L(r), and the set of left-invariant (resp., right- 
invariant) vector fields on G. If A is a section of t : AG —5- M, the corresponding 
left-invariant (resp., right- invariant) vector field on G will be denoted X~ (resp., X), 
where 

= (T e{m) l g )(X(j3(g))), (2.1) 
- -(T e(a(s)) r s )((r e(a(s)) i)(X(a( 5 )))), (2.2) 

for jeG. 
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Using the above facts, we may introduce a Lie algcbroid structure ([•, on 
AG, which is defined by 

{x,Y] = [1c,P], p(X)(x) = (T e{x) 0)(X(x)), (2.3) 

for X, Y G r(r) and x G M. We recall that a Lie algebroid A over a manifold 
M is a real vector bundle r : A — > M together with a Lie bracket [•, •] on the 
space r(r) of the global cross sections of r : A — > M and a bundle map, called the 
anchor map, such that if we also denote by p :T(t) —¥ X(M) the homomorphism 
of C°° (M)-modules induced by the anchor map then 

lX,fY\ = fiX,Y\ + p{X){f)Y, 

for X,Y £ T(r) and / G C°°(M). The triple (A, [•, -],p) is called a Lie algebroid 
over M (see [13]). 

Next, we will present some examples of Lie groupoids and their associated Lie 
algebroids which will be useful for our purposes. 

Any Lie group G is a Lie groupoid over {e}, the identity element of G, and the 
Lie algebroid associated with G is just the Lie algebra g of G. On the other hand, 
given a manifold M, the product manifold M x M is a Lie groupoid over M, called 
the pair or banal groupoid, in the following way: a is the projection onto the 
first factor and /3 is the projection onto the second factor, m((x, y), (y, z)) = (x, z), 
for (x,y), (y, z) G [M x M) 2 , e(x) = (x,x), for all x £ M, and i(x,y) = {y,x). 
The Lie algcbroid A(M x M) of this groupoid is isomorphic to the tangent bundle 
t m : TM -> M. 

The cotangent groupoid. (See [5], for details). Let G =4 M be a Lie groupoid. 
If A*G is the dual bundle to AG then the cotangent bundle T*G is a Lie groupoid 
over j4*G. The projections (3 and <S, the partial multiplication ©t*Gi the identity 
section e and the inversion I are defined as follows, 

P(p 9 )(X) = p g ((T tWg)) l g )(X)), for p g G T g *G and X G ^ (ff) G, 
«K)(T) = ^((T e(Q( , l)) n l )(T-(T £(Q(/l)) (6o/3))(y))), 

for G X£G and Y G A a ( h )G, 
(p g ©t-g f/ l )(T (s!fe) m(X s , n)) = p g (X g ) + ^(Fh), 

for (X 9! n) €T (gth) G 2 , 
e{p x ){X t(x) ) = p x {X t(x) - (T e{x) (e a a))(X e ^))), 

for p x G A* X G and G T e{x) G, 

l{p g ){X g -,) = -u g {{T g -n){X g -i)), for p g G T ff *G and Jf fl _, G T g -iG. 

Note that e(^4*G) is just the conormal bundle of M = e(M) as a submanifold of 
G. In addition, A*G is endowed with a linear Poisson structure. For this Poisson 
structure on A* G, a : T*G -> A* G is a Poisson map and /3 : T*G — > yl*G is an 
anti-Poisson (for more details, see [51 IT5] ^) . 
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3. Discrete dynamics in implicit form 

Motivated by the studies of implicit differential equations in P31 [T5J IS] and its 
generalization to Lie algebroids in , we introduce the notion of implicit difference 
equations on Lie groupoids and the corresponding notion of integrability. 

Definition 3.1. An implicit difference equation on a Lie groupoid G =} M is 
a submanifold E of G. 

An admissible sequence on the Lie groupoid G =4 M is a mapping 7g : I n 
Z — > G such that (3(jQ(i)) = a(7(j(i + 1)) for all i,i + 1 G / HZ. Here L is an 
interval on R. 

A solution of an implicit difference equation E C G is an admissible sequence 
7g : / fl Z — > G on G such that 7g(*) G E, for all i G I HZ. 

Let £ be a submanifold of the Lie groupoid G =4 M. Then E is said to be 

1) forward integrable at g G E if there is a solution 7g : Z + — > E C G 
with 7G (0) = .9. Here, Z+ = {n G Z | n > 0}. 

2) backward integrable at g G E if there is a solution 7c : Z~ — > E C G 
with 7G (0) = .g. Here, Z" = {n G Z | n < 0}. 

3) integrable at g E E if if there is a solution 7g : Z — > E C G with 
7g(0) = ff, that is, it is backward and forward integrable at g G E. 

If these conditions hold for all g, we say that E is forward integrable, backward 
integrable or integrable, respectively. 

Proposition 3.2. Let E be a submanifold of the Lie groupoid G =4 M. Then, 

1) E is forward integrable if and only if for each g G E exists at least an h G E 
such that (g,h) G G2 or, equivalently, E C f3~ 1 (a(E)). 

2) E 1 backward integrable if and only if for each g G E exists at least an hi G E 
such that (h',g) G G2 or, equivalently, E C a~ 1 (/?(£')). 

3) £ is integrable if and only if E C a" 1 (/?(£)) n /^(ali;)). 

Proof. I) Suppose that _E is forward integrable. Then, for all g G -E there is an 
admissible sequence 7g : Z+ — > E such that 7g(0) = g. So, we have that h = 
7g(1) G E satisfies (g,h) G G 2 . 

Conversely, if g G -E there exists & g\ € E such that (g,gi) G G2. Since 51 G E 
and using again the hypothesis, there is (72 G E satisfying (51,92) G G2. Therefore, 
we can construct a sequence {i7i}igz+ C E where go = g and (s'i , ) G G2. As a 
consequence, the sequence 7g : Z + ->£, m 7g(*) = ffij is a solution of E with 
7g(0) = fl. 

With a similar argument we deduce 2), and 3) is a direct consequence of 1) and 
2). ■ 

Example 3.3. Let M be a manifold and ip G C°°(M) be a diffcomorphism on M. 
Consider the implicit difference equation E v of the pair groupoid M x M given by 

E v = {(x, ip(x)) G M x M I x G M}. 
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Then, if (xq, ip(xo)) S E v , the map 7 : Z — > E v dchncd by 

1 (t)^(^(x ),^ +1 (x )), forieZ, 

where (p~* = (ip^ 1 ) 1 , is a solution of the implicit difference equation E v . Therefore, 
E v is integrable. 

This example can be generalized to bisections on Lie groupoids. We recall that 
a bisection on a Lie groupoid is a map o : M — » G which is a section of a and such 
that j3 ocr is a diffeomorphism. Given a bisection, a, if = <r(M) is the associated 
implicit difference equation on G then it is integrable, since for any g = <j{x) G E ai 
the sequence 

1 (i) = a{{Poa) l {x)) 1 ioiieZ, 
is a solution with 7(0) = g. 

The sets 

Ef = {g E E \ E is forward integrable at g}, 

Eb = {9 £ E I E is backward integrable at 5}, 

= {9 ^ E \ E is integrable at g}, 

are called the forward integrable, backward integrable and integrable parts 

of E, respectively. The implicit difference equation E is integrable (rcsp. forward 
integrable or backward integrable) if E = Efb (resp. E = Ef or E = Eb). Note 
that if E C E' then E f C E' f , E b C E' b and £/ b C E 1 ^. 

Proposition 3.4. Let E be a submanifold of a Lie groupoid G =t M. Assume that 
Efb, the integrable part of E, is a submanifold. Then, Efb is integrable. 

Proof. If g G Efb then there exists an admissible sequence 70 : Z — > E C G with 
7g(0) = 3- We have that 7g(Z) C Efb- Indeed, for all n £ Z there is an admissible 
sequence ^c.n ■ Z — > B given by 

lG,h{i) — -f G (i + h), for all i G Z. 



A similar result holds for Ef and Eb- Thus, Ef and Eb are forward and backward 
integrable, respectively. 

Now, it is easy to prove the following 
Proposition 3.5. Let E be a submanifold of a Lie groupoid G =t M. 

(i) If Ef is the forward integrable part of E and E' is a submanifold of G such 
that Ef C E' C E, then Ef is the forward integrable part of E' . 

(ii) If Eb is the backward integrable part of E and E' is a submanifold of G 
such that Eb C E' C E, then Eb is the backward integrable part of E . 

(iii) If Efb is the integrable part of E and E' is a submanifold of G such that 
Efb C E' C E, then Efb is the integrable part of E' . 
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From Proposition ^. 21 we have the following obvious relations 
E f C ^WEf^Q^HE)) 
E b C a-^PiE^Ca-^PiE)) 
E fb C /3- 1 (a(^ /6 ))na- 1 (/3(.B /6 )) C/3- 1 (a( J B))na- 1 (/3(.B)) 
Therefore, 

C Enf'fafB)) 

^ C £na^(^( £ )) 
E fb C Bnj3 _1 (a(B)) n a" 1 (/?(£)) 

Propositions 13.21 and 13.51 suggest a method for extracting the integrable part of 
an implicit difference equation under some regularity hypothesis. 

Let E C G be an implicit difference equation of a Lie groupoid G =4 M. 

For forward integrable part. Consider the sequence of sets of M 

C° f = a(E), C) = a(E n /T^G?)), . . . , C) = a(E n /T^G*" 1 )), • • • 

and the sequence subsets of G 

E° f =E, E}=En ^(Cj), ...,E)=EC\ /T^G)" 1 ), . . . 

Observe that a(E k f ) = G), C) C G^ 1 and C E^" 1 . 

For backward integrable part. Consider the sequence of sets of M 

D° b = P(E), D\ = (3{E n a~\D° b )), ...,D%=p(En a" 1 ^" 1 )), . . . 

and the sequence of subsets of G 

E° b =E, El = EH a-^Dl), ... 7 E k =E(l a" 1 ^ 1 ), . . . 

Also then D$ C D^ 1 and E\ C E^ 1 . 

For integrable part. Consider the sequence of sets of M 

C° fb = a(E), D% = p(E), 

c) b = a(Ena- 1 (D fb )nr 1 (C°f b )), D^pn^^Jnr 1 ^)) 

c% = ^na-^^jnr 1 ^ 1 )), D^^Bna-^nr 1 ^- 1 ))- 

and the sequence of subsets of G 

£° 6 =E,E} b = EH a-HD%) n r 1 ^), • . . 
E) h = En a-\D)^)) n p-HC^ 1 ), . . . 

Then C% C Cjj" 1 , £>) 6 C D^ 1 and E* 6 C 

It may happen that after a finite number of steps, the sets in two consecutive 
steps in the different sequences are equal: 

Cf = C k / +1 = C f , D k b » = D k b b+1 = D b , 

C fb ~ C fb V fb - V fb =JJ fb- 
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Then 

It is straightforward to see that if this happens then the subsequent steps in the 
different sequences are all equal. 

Moreover, observe that a(Ef) = Of, (5(E b ) = D b , a(Ef b ) = Ofb and ft (Efb) = 
D fb . 

Theorem 3.6. The forward integrable part of E is Ef, the backward integrable part 
is E b and the integrable part of E is Efb- 

Proof. Let g £ E f = E k / = E k / +1 . Then, since E k / +1 = E n (5~ 1 (C k /) 1 we have 

that (5(g) £ CJ = a(Ef), that is, there exists h G E kf+1 = E kf = Ef such that 
(5(g) = a(h). Therefore, using Proposition ^. 21 we have that E is forward integrable 
at g. 

Conversely, suppose that E is forward integrable at g. Then, there exists an 
admissible sequence 70 : Z + -> £ C G such that 7g(0) = g. Denote 7g(*) by gi, for 
all ieZ+. Then, (5(g kf ) = a(g kf+ i) and g kf £ E n (5~ l (C Q f ) = Ej, which implies 
that a(g kf ) € Cj, Now, since (5(g kf -i) = a(g kf ), we deduce that gk f -i & Ej. If we 
repeat this argument kf times, we conclude that g = g £ E^ = Ej = Ef. 

The corresponding results for backward intcgrability and intcgrability are proved 
in a similar way. ■ 

Remark 3.7. We have to stress that the concept of implicit difference equation 
can be to generalized to the setting of a pair of manifolds G and M endowed with 
two surjective submersions a, (5 : G — > M . The corresponding notions of admissible 
sequences, solutions, integrable equations and the integrability algorithm can be 
straightforwardly generalized. This more general approach will be useful in Section 

mi o 

Remark 3.8. For some concrete applications (see, for instance, the next example) 
instead of a unique subset E of the Lie groupoid G =4 M, we have a sequence of 
subsets {£'fc}o<fe<oo- A generalization of this type is also useful for discretization 
of time-dependent systems. It is not hard to extend our setting for this situation 
and also to extract the integrable part of a sequence of this type. For example, the 
forward integrable part is extracted using the following procedure. Consider the 
/c-sequence of sets of M defined for k > 0: 

[C k ]° f = a(E k ), [C k ]) - a(E k n(5~ 1 ([C k+ i] a f )), [C k ) f = a^n/T 1 ^([C fe+1 ]}- 1 )), 
and the subsets 

[E k ] f = E k , [E k ]\ = E k n (5-\[C k+ y]% [E k ]\ = E k n ^([Ckff 1 ). 



If the algorithm stabilizes in a sequence of subsets {[i^fc]^} we will call this set the 
backward integrable part of the sequence {E k }. o 
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Example 3.9. Finite difference methods for time-dependent linear diffe- 
rential-algebraic equations The proposed algorithm is also useful for numerical 
methods for differential-algebraic equations (DAEs). These systems of equations 
appear in a great number of applications: mechanical systems with constraints, 
electrical networks, chemical reacting systems... See, for instance, [3J [52]. These 
types of systems take the form F(t,x,x) = including as a particular example 
the explicit differential equations. The system that we will study are the so-called 
time-dependent linear differential-algebraic equations: 

A(t)x + B(t)x = b(t), A(t),B(t) g£(R™,R"), b(t) G R n , t G R 

Consider an explicit Euler approximation to this system 

x k+1 -x k 

A k 7 1- B k x k = b k 

h 

where h > is the step size and A k = A(t +kh), B k = B(t +kh) and b k = b(to+kh) 
and xo,to initial conditions. Define the subsets of the pair groupoid R™ x R n 

E k = {(x k ,x k+1 ) £ M" x 1" | A k Xk+ \~ Xk + B k x k = b k } 

h 

We consider V k the image of the linear map associated with A k . Then, V k = R™ (in 
this case, x k+ \ = x k + hA^ 1 ^ — B k x k )) or there exists a matrix Q k G £(R n ,R") 
such that the kernel of its linear map is just V k , i.e., Q k A k = 0. Then, one may 
easily prove that 

[C k ]° f = {xeR n | Q k B k x - Q k b k = 0} . 

Therefore, 

[E k ]) = {(x k ,x k+1 ) e R" x R" | A k Xk+1 ~ Xk +B k x k = b k , 

Qk+\B k+ ix k+ i — Q k +\b k+ i = 0} 
= {(x k ,x k+1 ) G R" x R" | Q k+1 B k+1 x k+1 - Q k+1 b k+1 = 0, 
(A k + Q k+1 B k+1 )x k+ i = {A k - hB k )x k + hb k + Qk+ibk+x, } ■ 

It is possible to continue the algorithm, considering a([i?fc]y) = [C^Jj, but in many 
cases of interest (as, for instance, discretizations of index one DAEs) we have that 
the matrix A k + Q k+ \B k+ \ is regular, then 

x fe+ i = (A k + Qfc+i-B/c+i) -1 {{A k - hB k )x k + hb k + Q k+ ib k+1 ) . 

In this case, a([E k ]j) = [C k ]j = [C k ]^ and the algorithm stops in this step (see [53] 
for more details). 

4. Implicit systems defined by discrete Lagrangian functions 

4.1. Implicit systems defined by discrete Lagrangian functions. Lagran- 
gian submanifolds. Let L : G — ► R be a discrete Lagrangian on a Lie groupoid 
G =t M . Then, Sl = dL{G) C T*G is clearly a Lagrangian submanifold of the sym- 
plectic manifold T*G and Sl is also an implicit difference equation on the cotangent 
groupoid T*G =4 A*G. If we apply the forward integrability algorithm to Sl, then 
there is an interesting relation with the discrete Euler-Lagrange equations for L (see 
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(T51 |2"5] and the Appendix [A] for more details on discrete Lagrangian Mechanics on 
Lie groupoids) , which we describe in this section. 

Proposition 4.1. Let G =4 M be a Lie groupoid and L : G — > R be a discrete 
Lagrangian on G. The first step, (Sl) 1 /, to obtain the forward integrable part of 
Sl = dL(G) are the points dL(g) G Sl such that there exists an element h G G, 
satisfying (g,h) G G*2 and, in addition, (g,h) is a solution of the discrete Euler- 
Lagrange equations, that is, 

d[Lol g + Lor h oi] (e(x))\ AxG = 0, 

where x = (3(g) — a(h). 

As a consequence, if (resp. a) is the target (resp. source) map of the cotangent 
groupoid T*G =4 A*G defined in (|2.4[) , then 

(Sl)/ = {is-, h) 6 Gi | (g, h) is a solution of the discrete Euler- Lagrange eqns.} 
= {(g,h)£G 2 \p(dL(g))=a(dL(h))} 

Proof. Using the definition of the source map a : T*G -> A*G (see (|2.4|) ) and the 
fact that 

X - T e{x) (e o 0)(X) = -T e{x) i(X), for X G A X G, 

we have that 

C° f = &(S L ) = {a* G A*G\a* = -d{Lor h oi)(e{a(h))) G A* a(h) G, for some h G G}. 

Therefore, (Sl)j (/?) _1 (C°) is the set of points dL(g) G Sl such that there 

exists h G G, with (g, h) G Gi and 

d [L o l g + L o r h o i] (e(i))u iG = 0. 



Let L : G — > R be a hyperregular discrete Lagrangian function (see Appendix lAf. 
This implies that the discrete Lagrangian evolution operator T l ■ G — > G is given by 
T l = (¥-L)- 1 o (W+L). Then, if g € G, the map T T * G : Z ^ S L = dL(G) C T*G 
defined by 

T T . G (i) = dL(Ti(fl)), forteZ, 
is a solution of the implicit difference equation Sl and Tt*g(0) = dL(g). 
Thus, we have proved the following result. 

Corollary 4.2. Let G =3 M be a Lie groupoid and L : G — > R be a discrete 
hyperregular Lagrangian function on G. Then, the implicit difference equation Sl = 
dL(G) is integrable. 

Under the same hypotheses as in Corollary 14.21 we may define the map 

$ : G G x G (F ^T L) A*G x A*G 

g^(g,T L (g)) .— > (¥-L(g),¥+L(g)) 

A direct computation, using (|A.3|) . shows that $(G) = (a,/3)(5i). 
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On the other hand, if f L = F+L o (F~L) -1 : A*G —¥ A*G is the discrete 
Hamiltonian evolution operator then one observes that 

graphf L = $(G) = (aJ)(S L ). 

Using that Sl is a Lagrangian submanifold of T*G and that (<5,/3) : T*G — > 
A*G x A*G is a Poisson map, A*G denoting A*G endowed with the linear Poisson 
structure changed of sign, then (a, (3) (Sl) is a coisotropic submanifold of A* G x A* G 
(this is a particular case of [55J Corollary 2.2.5]). 

From the previous discussion, we can conclude the following result (see also [29]). 

Proposition 4.3. Let G ^ M be a Lie groupoid and L : G — >• R be a discrete 
hyperregular Lagrangian on G. Then, the discrete Hamiltonian evolution operation 
T l preserves the Poisson structure on A*G. 

4.2. Example: the pair groupoid. Let Q be a manifold and L : Q x Q —> R 

be a discrete Lagrangian. We know that dL(Q x Q) is a Lagrangian submanifold 
of T*(<2 x Q). On the other hand, T*(Q xQ) 4 T*Q is a symplectic groupoid 
with the canonical symplectic structure of T*(Q x Q). The structural maps of 
T*(Q x Q) =X T*Q are 

P:T*(QxQ)^T*Q, ( 7go , 7gi )^ 7gi , 
a:T*(QxQ)^T*Q, ( lqo , lqi ) -> - 7go , 

m : T*(Q x Q) 2 -> T*Q, ((790.791), (-7?i,7g 3 )) ^ (790,792), 
i : T*(Q xQ)-> T*(Q x Q), ( lqo , lqi ) h> (- 7(?1 , - 7 , ), 
£:T*Q^r(QxQ), 79^(79,-79)- 

Moreover, 

(3(dL(q ,qi)) = D 2 L(q ,q 1 ), 
a(dL(q 1 ,q 2 )) = -D 1 L(q 1 ,q 2 ), 

for (90, 9i), (91, 92) € Q x Q. Thus, the discrete Eulcr-Lagrange equations for the 
points (qo, 9i)j (91, 92), that is D 2 L(qo, 91) + D\L(q\, 92) = 0, are equivalent to 

P(dL(q ,q 1 )) = a(dL(q 1 ,q 2 )). 

Therefore, we obtain the following result (which can also be seen as a consequence 
of Proposition 14. ip 

Corollary 4.4. Let Q be a smooth manifold and L : Q x Q — > R be a discrete La- 
grangian function. Then, there exists a bijective correspondence between composable 
points in dL(Q x Q) and solutions of the discrete Euler- Lagrange equations. 

Remark 4.5. T*(Q x Q) =4 T*Q is a Lie groupoid which is symplectomorphic 
to the pair groupoid T*Q x T*Q =4 T*Q when we consider the symplectic form 
(-Oq,Oq) on T*Q x T*Q. 
The isomorphism is given by 

*:T*(QxQ) T*Q x T*Q 

(790, 791) ^ (-790: 7 9 i) 
Note that o dL)(q , 91) = {-D 1 L{q ,q 1 ),D 2 L{q n ,q 1 )). o 
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Now, suppose that L : Q x Q —> R is hypcrrcgular. Then 

T L = (W-L)- 1 o¥+L : Q x Q -> Q x Q (4.2) 

is a discrete Lagrangian evolution operator. In addition, we have that the image of 
the map 

QxQ {M 4 L) (QxQ)x(Qx Q) rL ^ L T *Q x T*Q 
(<7o,<7i) !-> ((Qo,Qi),"^l(qo,Qi)) ^ (-D 1 L(q ,q 1 ),D 2 L(q a ,q 1 )) 

is just the Lagrangian submanifold o dL(Q x Q). 

On the other hand, we can consider the discrete Hamiltonian evolution operator 
T l = V+L o (F-L)- 1 : T*Q -> T*Q. A direct computation proves that 

*(Q x Q) = (* dL)(Q x q) = (7d, f l)(T*Q). 

We remark that, since : T*Q — > T*Q is a symplcctomorphism, we deduce 
(Id,f L )(T*Q) is a Lagrangian submanifold of (T*Q x T*Q, (-Qq, Oq)). 

4.3. Example: a discrete singular Lagrangian. Consider the discretization of 
the degenerate Lagrangian L : TR 2 — > R given by L(cc, y, i, y) = \x 2 + \x 2 y (see 
[2]). Take, for instance, a consistent discretization of this lagrangian: 

Lj:R 2 xl 2 ^l, Lk( Xljl/lia . 2jl/a ) = I^LZ£i^ +ia; 2 2/1 

Then S l h C T*(R 2 x R 2 ) defined by 

c X2 - xi 1 2 a; 2 -£i m , , > to4 

6 L h = <^ (asi, y!,X2,V2, ^ h xij/i, — p — >°) I (^1,2/1, £2, 2/2) eR 

The algorithm produces the sequence: 

C) = 5(S lS )=fe^-^-^)l(i,!/,i)eK 3 )}cn 2 



Xl I ™ „. I I n\ „. \ r- m>3 



S'l^J^. = 0,y 2 ; -p- + xiyi, -^x, -^2" , 0) | (x 1 ,yi,y 2 ) G 

C 2 = a((5 L5 )p = {(x,y;-^-^,-ix 2 )|(.T, 2/ )eR 2 )}cT*R 2 



(%) = {(0, yi ,0,j/ 2 ;0,0,0,0) I 6 R 2 } 



9 

Cf = 5((5 LS ) 2 ) = {(0, 2 /;0,0)|yGR)}cT*R 2 



Since f<$xj»^ = (^L^j then the algorithm terminates with (s L h\ = f$z,J>J > 
which is the forward integrable part of the implicit difference equation S L u . 
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4.4. Implicit systems defined by discrete nonholonomic Lagrangian sys- 
tems. Let (L, M c , D c ) be a discrete nonholonomic Lagrangian system. Consider 
the submanifold S(L,M a ) T*G given by 5'(£ i m c ) = dL(M c ), that is, 

S(L,Ma) = { dL (*M c (g)) g £ M c } . 

Composing the source and the target maps a and j3 of the cotangent groupoid 
T*G =4 A*G as defined in Equation (JH) with i*^ : AG* -> T>*, the dual map of 
the canonical inclusion i^ c : D c — >• AG, we have the maps 

5a,. = i$, o o a : T*G -> D; 

0d c =ii o o^:T*G^D*. 1 J 

If we apply the forward integrability algorithm to the implicit difference equation 
S(l.m c ) m T*G, but changing the maps a and /3 for the new maps an c and /3u c 
(see Remark l3.7p . then the first step of the algorithm yields a close relation with the 
discrete non-holonomic Euler-Lagrange equations. Indeed, if we have the following 
result. 

Proposition 4.6. Let G =t M be a Lie groupoid and (L,M C , D c ) be a discrete 
nonholonomic Lagrangian system on G. The first step, (S'(l ! m c ))/; t° obtain the 
forward integrable part of SfL t M.\ = dL(M c ) for the maps ai> c : T*G — > D* and 
Pd c : T*G — > D* c , defined in (|4.3[) . are the points dL(g) € S(l,m c ) such that there 
exists an element h £ M C; satisfying (g,h) <= G2 and, in addition, (g, h) is a solution 
of the discrete non-holonomic Euler-Lagrange equations, that is, 

d[Lol g + Lor h oi] {e{x))\(<r> c ) x = 0, 

where x — /3(g) — a(h). 
Ln other words, 

(S L )} = {(g,h)eG 2 n(M c xM c )\ 

(g, h) is a solution of the discrete nonholonomic Euler-Lagrange eqns.} 
= {(g,h) e G 2 n (M c x M c )\fo c (dL(g)) = av e (dL(h))} 

Proof. The proof is analogous to that of Proposition 14. 1[ taking into account the 
version of the discrete nonholonomic Euler-Lagrange equations given by Equation 
Ml . ■ 

4.5. Example: The discrete Chaplygin sleigh [71 |S]. To illustrate the results 
contained in the previous section, we will describe a discretization of the Chaplygin 
sleigh system (previously considered in [TUj). 

The Chaplygin sleigh system describes the motion of a rigid body sliding on a 
horizontal plane. The body is supported at three points, two of which slide freely 
without friction while the third is a knife edge, a constraint that allows no motion 
orthogonal to this edge (see HJ). 

The configuration space of this system is the group SE(2) of Euclidean motions 
of IR 2 . An element A £ SE(2) is represented by a matrix 




with 9,x,y £ E. 
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Thus, (9,x,y) arc local coordinates on SE(2). 

A basis of the Lie algebra se(2) = R 3 of SE(2) is given by 






ei = U U U , e 2 



and we have that 

[e,ei]=e 2 , [e,e 2 ]=-ei, [ei,e 2 ]=0. 
An element £ £ se(2) is of the form 

£ = w e + ui ei + v 2 e 2 
The exponential map exp : se(2) = R 3 — > SE{2) of SE(2) is given by 

sin cj cos u> — 1 . cos w — 1 . sin w . 

exp(w,ui,i> 2 ) = [0J,Vi \-v 2 ( ) + u 2 ), ifw^O, 

w w ww 

and exp(0, v±, w 2 ) = (0,?;i,w 2 ). The restriction of the exponential map to the open 
subset U =} — 7r,7r[x]R 2 C R 3 = se(2) is a diffeomorphism onto the open subset 
exp(J7) of SE(2). 

A discretization of the Chaplygin sleigh may be constructed as follows (see [TU] 
for more details): 

First of all, the discrete Lagrangian L : SE(2) — > R is given by 
L(A) = ^Tr (ASA T ) - Tr (A3), 

where J is the matrix: 

(J/2) + ma mab ma 

mab (J/2)+mb 2 mb 

ma mb m 

(see [S]). In terms of the coordinates (9, x, y), the discrete Lagrangian can be written 

as 

L(6, x, y) = (max + mby — ma 2 — mb 2 — J) cos 9 + m(ay — bx) sin 9 
+ ^(( x -a) 2 + (y-b) 2 )+\(J-m) 

Second, the vector subspace D c of the Lie algebra sc(2) is given by 

D c = span {e,ei} = { (uj,v 1 ,v 2 ) <E se(2) | v 2 = } . 

Finally, the constraint submanifold M c of SE(2) is M c = exp(C/ n D c ), that is, 

M c = { (6,x,y) e SE(2) \ -n < 6 < tt, 6 + 0,(1 -cos6»)a; - j/sin6» = 0} 
U{(0,x,0) G SE(2) 
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Let us calculate the set S^Mc) = dL(M, c ). We have that (9,x,y;pe,p x ,p y ) £ 
T*SE(2) = M 6 belongs to SfL,M c ) if it verifies the following equations: 

pe = m(ay — bx) cos 9 + (ma 2 + mb 2 + J — max — mby) sin 6, 

p x — ma cos 6 — mb sin 9 + m(x — a), 

p y — mb cos 9 + masm.9 + m(y — 6), 

with (9,x,y) <E M c . 

Moreover, using (|2.4[) . we deduce that the maps f) : T*SE(2) -> se(2)* and 
a : T*SE(2) — > se(2)* are just the pullbacks by the left and right translations, 
respectively. Now, the mappings 5<d c : T*G ->■ D* and /3d c : T*G -> CD* are: 

5d c (p e d6> + P^rfa; + p y <iy) = (pe - + ^)e* + p x e\ 

Pd c (ped9 + p x dx + p y dy) = p g e* + (p x cos 9 + p y sin 9)e\ 

where pod9 + p x dx + p y dy € x y ^G and {e*, e^} is the dual basis of {e, ei} and 
therefore span {e*, e\} = D*. 

Thus, the first step of the algorithm allows yields the set of points (9k,Xk,yk) £ 
M c , with k e {1, 2}, such that 

P>D c (dL(6i,xi,yi)) = aT> c (dL(9 2 ,x 2 , j/2)), 

which are just the discrete Euler-Poincare-Suslov equations, 

-ma cos 9\ — mb sin 9\ 4~ ma \ / mx 2 + ma cos 9 2 
Ymx\ cos 0i + m?/i sin^i y I —mbs'm9 2 —ma 

m(ayi — bx\) cos9\ — m{ax\ + byi) sin^i \ / may 2 ~mbx 2 
+ (ma 2 + mb 2 + J) sin #1 y I +(ma 2 + m6 2 + J) sin0 2 

where (9k,Xk,yk) & 3Vt c , with fc € {1,2}. 

4.6. An approach to implicit discrete Hamiltonian systems. Along the pa- 
per, we have focused our attention to the case of discrete lagrangian mechanics 
including this theory in the setting of discrete implicit systems. Of course, we can 
also adopt a dual point of view, that is, the hamiltonian formalism. In the continu- 
ous setting, the hamiltonian formalism is specified given a hamiltonian function H 
on A*G equipped with its canonical linear Poisson bracket. The dynamics is given 
by the associated hamiltonian vector field Xh- In the sequel, we will show how to 
obtain a lagrangian submanifold of the symplectic groupoid T*G =t A*G directly 
from the flow of Xh- 

Assume that G =4 M is a Lie groupoid with corresponding algebroid r : AG —> 
M. Then, we have already explained (see (|2.4[1 ) that the cotangent bundle T*G is 
a groupoid over A*G. Moreover, there is a Poisson structure on A*G (induced from 
the Lie algebroid structure on AC?) such that the source map a : T*G — > A*G (resp. 
the target map : T*G — > A*G) is a Poisson (resp. anti-Poisson) map. 

Let H : A*G — > R be a Hamiltonian function. Then, since A*G is a Poisson 
manifold, we have the hamiltonian vector field Xh- On the other hand, consider the 
pull-back of the 1-form dH £ £l 1 (A*G) by the source map a. Using the canonical 
symplectic structure on T*G, we have the corresponding Hamiltonian vector field 
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Xa*H- It is easy to see that both vector fields are 5-relatcd. Indeed, if Ha*g (resp. 
ILg) denotes the Poisson structure on A*G (resp. on T*G), since 5*(IIg) = ILa*g, 

X h = H A , G (dH) = (5,n G )«(dff) = a*((IL G )t(da*H)) = a*(X & * H ), 

where H^(dF) denotes contraction by dF, i.e., IT" (dF) (7) = H(dF, 7) for any 1-form 
7. Thus, we have that the following diagram is commutative 



T*G 



T*G 



a 



A*G 



H 



A*G 



where <p H (resp., tjj H ) is the flow of X^-h (resp., Xh)- In addition, since <p H : 
T*G — > T*G is a symplectomorphism, i : T*G — > T*G is an anti-symplectomorphism 
and e : A*G — > T*G is a Lagrangian embedding, we deduce that the submanifold 
S H defined by 



S H = {u g e T*G I u g = 
is a Lagrangian submanifold of T*G. 



i(tp H (e(x))),Vx g A*G) c T*G, 



Remark 4.7. We must note that the previous discussion is a particular instance of 
the construction of Lagrangian bisections on symplectic Lie groupoids (see [5]). o 



Example 4.8. Let H : T*Q 



be a Hamiltonian function on T*Q. In this case, 



using canonical coordinates (q,p), the Poisson structure on T*Q is IIq = — ^ A 
and the Hamiltonian vector field is given by 

dH d dH d 
H dp dq dq dp 

Now, let G = Q x Q be the pair groupoid. Then, on T*(Q x Q) the source map 



a : T*(Q x Q) 

Hqxq = 4A^ + ^ 7 A s f r ,we have 



T*Q is just (qo,qi,Po,Pi 
dH 



da*H 

Xa'H 



dq 



(qo,—po)- As a consequence, since 
dH , 



dqo 
dH d 

dp dq 

Finally, since e : A*G —> T*G, (q,p) >-} (q,q, —p,p), we have that 

{(q, (iPi) H (q,p), -p, (ih)* H (q,p)) \ (q,p) 6 T*Q} , 



9p 

ai? d 

dq dp 



where 



(q,p) H- ip H (q,p) = (W>i ) H (q,p), (tp2y H (q, P)j is the flow of Xjj. 



Remark 4.9. Let L : TQ ->Rbea regular continuous Lagrangian on the manifold 
Q. The exact discrete Lagrangian for a small time-step h > ([T7]) is the 
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given by 
h 

L(a i(t),a 01 (t))dt 



discrete Lagrangian L d : Q x Q — 



where (Joi(i) is the unique solution of the Euler-Lagrange equations for L which 
satisfies the boundary conditions coi(O) = qo and croi(h) = q\. 

The Legendre transformations, ¥^ h and YL, for Lh and L respectively, are related 

as follows (see [17] ) 

Y+L%(q , qi ) = YL(a 01 (h),& i(h)) 
Y-L%(q , qi ) = FL(a O i(0),doi(0)) 

Therefore, if L is (hyper-)regular then L h d is (hyper-) regular. 

Suppose that L is hyperregular. Denote by H : T*Q — > K the corresponding 
Hamiltonian function, i.e. 

// = BioFr 1 

El being the Lagrangian energy associated with L. 

In [17] it is proved that for a small time-step h, the following diagrams relate the 
flow ip^ : T*Q -x T*Q at time h of the Hamiltonian vector field X H G X{T*Q) and 
the discrete Lagrangian evolution operator T L h : Q x Q — x Q x Q associated with 
L h d defined as in (|4~2l 



QxQ 



QxQ 



QxQ 



F+ 



T*Q 



T*Q 



T*Q 



A direct computation, using this result, shows that 

= (* o dL h d )(Q xQ) = graph ^ = 



Here * : T*(Q x Q) -> T*(Q x Q) is the isomorphism described in (|4.1|) . Thus, 

It would be interesting to extend this result for a general regular Lagrangian 
L : AG —x K on the associated Lie algebroid ^4G of a Lie groupoid G =4 M. 
Previously, its is necessary to describe the exact discrete Lagrangian L d : G — x M 
induced by L (see Q15] for more details). 

The results of this section motivate that we may introduce, in the general context 
of a Lie groupoid G, the notion of an implicit discrete Hamiltonian system as a 
Lagrangian submanifold of T*G. 

o 
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5. Conclusion 

We have developed an algorithm which permits to extract the integrable part of 
an implicit difference equation. This algorithm produces a sequence of submanifolds 
(or generally subsets) which encodes where there exist well defined solutions of the 
discrete dynamics. As an application, our method allows us to easily analyze the 
case of discrete (nonholonomic) Lagrangian systems from an implicit point of view, 
in particular the situation of singular lagrangians, and we will expect that it will 
be an useful tool for discrete optimal control theory. We will discuss this topic in a 
future paper. 

Finally, in our opinion, the developments about implicit systems defined by 
Hamiltonian functions open doors to future research in the theory of generating 
functions for Poisson morphisms (using standard symplcctic techniques) , the theory 
of discrete hamiltonian systems in a Poisson context [1], the Dirac theory of con- 
straints in the discrete setting, a discrete Hamilton- Jacobi theory... among other 
research topics. 

Appendix A. Discrete Lagrangian Mechanics on Lie groupoids 

A discrete Lagrangian system consists of a Lie groupoid G =4 M (the dis- 
crete space) and a discrete Lagrangian L : G — >• ML 

For g £ G fixed, we consider the set of admissible sequences: 
Gg ={jg '■ [lj N] H Z — >• G | 7g is an admissible sequence and 7g(1)- -7g(A) = g } . 
We may identify the tangent space to at 7g with 

T 7G Cf = {(vi,V2, • ■ ■ ,v N -i) | v k G A Xk G and x k = ^ G (k)), 1 < k < N — 1} . 
(vi, i>2, • ■ • , vn—i) is called an infinitesimal variation of 7g. 

Now, we define the discrete action sum associated to the discrete Lagrangian 
L : G -> E by 

N 

SL( 7G )=^L( 7G (fc)). 
fc=i 

Hamilton's principle requires that this discrete action sum be stationary with respect 
to all the infinitesimal variations. This requirement gives the following alternative 
expressions for the discrete Euler- Lagrange equations (see |18|): 

1c(g H )(L) - Jt(g k+1 )(L) = 0, (A.l) 

for all sections X of r : AG —> M. Alternatively, we may rewrite the discrete 
Euler-Lagrange equations as 

d[Lol gk +Lor 3k+1 oi] (e(a*))u fcG = 0, 

where (3{g k ) = a(g k+1 ) = x k - 

Thus, we may define the discrete Euler-Lagrange operator -Ddel^ : G2 — > 
A*G from G 2 to A*G, the dual of AG. This operator is given by 

D D ELL(g, h) = d [L o l g + L o r h o i] (e(x))i 



DISCRETE DYNAMICS IN IMPLICIT FORM 



19 



with j3(g) — a(h) = x. 

Let T : G — > G be a smooth map such that: 

- graph(T) C G2, that is, (g, T(g)) G G2, for all g G G (T is a second order 
operator) . 

- (g, T(g)) is a solution of the discrete Euler-Lagrange equations, for all g G 
G, that is, (D DEL L)( 5 , T(g)) = 0, for all g € G. 

In such a case 

X r (g)(L)-X~(T(g))(L)=0, (A.2) 

for every section X of AG and every g € G. The map T : G — > G is called a discrete 
flow or a discrete Lagrangian evolution operator for L. 

Given a Lagrangian L : G — > R we define the discrete Legendre transforma- 
tions F~L : G -> A*G and F+L : G -» A*G by 

(F"L)(/i)(w £ ( Q (, l ))) = -v £ ( a ( h ))(L o r h o i) 

= a(dL(h))(v e ( a ( h)) ), for u e(a(/l )) £ A q(/i )G, 

(A.3) 

(F + L){g){v t{l3{g}} ) = Vj(f){ g) ){Lol g ) 

= P(dL(g))(v e(m) ), for v t{fi(g)) G A^G. 

A discrete Lagrangian i : G — >• M is said to be regular if and only if the Le- 
gendre transformation ¥ + L is a local diffcomorphism (equivalently, if and only 
if the Legendre transformation ¥~L is a local diffeomorphism) . In this case, if 
(<?o> G G x G is a solution of the discrete Euler-Lagrange equations for L then, 
one may prove (see [TS]) that there exist two open subsets Uq and Vb of G, with 
.go G C/o and /io G Vb, and there exists a (local) discrete Lagrangian evolution oper- 
ator Tl ' Uq — > Vq such that: 

(i) T L (.g ) = ho, 

(ii) Tl is a diffcomorphism and 

(iii) Tl is unique, that is, if Uq is an open subset of G, with g G f/gi an d 
T' L : Uq — > G is a (local) discrete Lagrangian evolution operator then 

"Filt/ont^ = ^z,|i7 n(7^- 

Moreover, if F + L and ¥~L are global diffcomorphisms (that is, L is hyperregular) 
then T L = (F-L)- 1 o F+L. 

If L : G — >• R is a hyperregular Lagrangian function, then pushing forward to A*G 
with the discrete Legendre transformations, we obtain the discrete Hamiltonian 
evolution operator, Tl '■ A*G — > A*G which is given by 

f L =F ± LoT l o(F ± L)- 1 = F+Lo (F^L)- 1 . (A.4) 

Appendix B. Discrete nonholonomic Lagrangian systems 

A discrete nonholonomic Lagrangian system on a Lie groupoid G =4 M 
is a Lagrangian discrete L : G — > R, a vector subbundle D c (the constraint 
distribution) of the Lie algebroid AG of G and a discrete constraint embedded 
submanifold M c of G such that dimM c = dimD c . 
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Let (X,M C ,2) C ) be a discrete nonholonomic Lagrangian system. The discrete 
nonholonomic Euler- Lagrange equations for the system (L d ,M. c ,D c ) are 
given by 

d(L d ol g + L d o r h o i)( e ( a; ))|(D c ) ai = 0, 
for (g, h) g G2 D (M c x M c ), with (3(g) = a(h) = x (for more details on this section, 
see [10]). 

For a discrete nonholonomic Lagrangian system, we can define the discrete non- 
holonomic Legendre transformations 

¥-(L,M c ,D c ) : M c -»• V* and ¥ + (L, M c , D c ) : M c — > D* 

as follows: 

¥^{L,M c ,D c )(h)(v e{a{h)) ) = -v £{a{h)) (Lor h oi), 

= a(dL(h))(v e ( a ( h ))), for e D c (a(/i)),(B.l) 

W + (L,M c ,T> c )(g)(v eWg)) ) = v eWg)) (Lol g ), 

= p(dL(g))(v e{m) ), for v e{m) £ D c (/3(g)). (B.2) 

If F _ L d : G -> AG* and F+L d : G -> AG* are the standard discrete Legendre 
transformations associated with the Lagrangian function L defined in Equation (|A.3[) 
and i? D : AG* — > D* is the dual map of the canonical inclusion i<r> c : D c —> AG 
then 

W-(L d ,M c ,'D c )=i*v c oW-L d oi Mc , ¥ + (L d , M c , D c ) = i* Dc o ¥ + L d o i Mc . (B.3) 

Using Equations (|B.1[) and (|B.2[) . the discrete nonholonomic Eulcr-Lagrange equa- 
tions are equivalent to 

¥-(L d ,M c ,D c )(h)=¥+(L d ,Mc,'r>c)(9)- (B.4) 
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